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According to the standard model, top quarks can be produced at the Tevatron pp collider 
via two mechanisms. One involves a virtual gluon that decays via the strong interaction into 
a ti pair. This mode has been observed by the CDF and D0 collaborations [|1],|2| . The second 
mechanism involves the electroweak production of a single top quark at a Wtb vertex, where 
W and b refer to the W boson and b quark. There are three ways of producing a single top 
quark: an s-channel process q'q—*tb, a t-channel mode q'g^tqb, and a final state generated 
via both the s- and t-channels, bg—>tW. For a top quark of mass 174.3 ± 5.1 GeV the 
predicted cross sections are 0.75 ± 0.12 pb for tb production || and 1.47 ± 0.22 pb for tqb 
production |5|], both calculated at next-to-leading-order precision with ^/s = 1.8 TeV. For 
tW production, the cross section is 0.15 pb at leading order ||. The errors on the cross 
sections include contributions from the choice of parton distribution functions (PDF), the 
uncertainty on the gluon distribution, the choice of scale, and the experimental error on the 
top-quark mass measurement. In this Letter, we use the notation u tb" to refer to both tb 
and the charge-conjugate process ib, and "tqb" to refer to both tqb and tqb. 

The D0 collaboration recently published its first results on the production of single top 
quarks 0. The search used a classical selection method to optimize the signal significance, 
based on the expected kinematic properties of the events. The 95% confidence level upper 
limits on the cross sections were determined to be 39 pb for the s-channel process and 58 pb 
for the t-channel process. The events contained an isolated electron or muon ||, missing 
transverse energy, and jets. At least one jet in each event was required to contain a "tagging" 
muon ||, used as an indication that the jet originated from the hadronization of a b quark. 
We have now developed a new and more powerful technique using arrays of neural networks 
that allow us to utilize the far more numerous untagged events in the search, as well as to 
improve the sensitivity for tagged events. This Letter describes the new method of event 
selection and presents significantly improved upper limits on the single-top-quark production 
cross sections. 



The D0 detector |TIJ in Run 1 (1992-1996) had three major components: a drift- 
chamber-based central tracking system that included a transition radiation detector, a 
uranium/liquid- argon calorimeter with a central module (CC) and two end calorimeters 
(EC), and an outer muon spectrometer. For the electron channel, we use 91.9 ± 4.1 pb -1 
of data collected with a trigger that required an electromagnetic (EM) energy cluster in 
the calorimeter, a jet, and missing transverse energy (Et)- For events passing the final 
selection criteria, the efficiency of this trigger is (90-99)%, depending on the location of 
the EM cluster in the calorimeter and on the presence of a tagging muon. In the muon 
channel, we use 88.0 ± 3.9 pb -1 of data acquired with several triggers that required either 
$Ti or a muon and a jet. The combined efficiency of these triggers is (92-98)%. A third 
data sample, obtained with a trigger requiring just three jets, is used for measuring one of 
the backgrounds. Since the multijet cross section is very large, this trigger was prescaled, 
and we have 0.8 pb -1 of such data. Each of the three samples contains approximately one 
million events. We reconstruct the events offline by applying the same criteria to identify 
electrons, jets, and isolated and tagging muons, as described in Ref. 0. 

Single-top-quark events have a readily identifiable final-state topology. The top quark 
decays to a W boson and a b quark. The W boson decays to a central (i.e., low pseudorapidity 
\r)\ ||11||), isolated electron or muon with high transverse energy and momentum (Et, Pt), and 
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a central, high-.&r neutrino. We infer the presence of the neutrino from the vector imbalance 
of Ex in the event. For an s-channel tb event, there are two central, high-E^ b jets, whereas 
for a t-channel tqb event, there is only one such b jet, plus a forward, light-quark jet, and 
a central \ow-E T b jet. About 11% of the time [fl^| , one of the b jets contains a tagging 
muon. Unfortunately, events like this are easily mimicked by many background processes. 
Before applying the initial selection criteria, the largest background in the electron channel 
comes from multijet events in which a jet is misidentified as an electron. We call this false 
isolated-electron background "mis-ID e." In the muon channel, the dominant component 
of the background is from multijet production with a coincident cosmic ray or beam-halo 
particle misidentified as an isolated muon. We call this false isolated-muon background 
"cosmic"; it is not included in the background model. Other backgrounds such as PU+jets, 
ti pairs, and bb pairs, contribute at the few percent level. The bb background affects only 
the muon channel, when a muon from a b decay is misidentified as isolated. We call this 
false-isolated muon background "mis-ID /^." 



TABLE I. Initial selection criteria. 

Pass triggers and online filters No mismeasured muons 

Exactly one good isolated lepton No mismeasured Et 

Any number of good tagging muons No mismeasured jets 
No photons 



Analysis Channel 



Variable 


e+jets/notag 


e+jets/tag 


/x+jets/notag 


//+jets/tag 


E T (e) 


> 20 GeV 


> 20 GeV 






\r,(e)\ 


< 1.1 (CC) 


< 1.1 (CC) 








> 1.5, < 2.5 (EC) 


> 1.5, < 2.5 (EC) 






pr(isol fi) 






> 20 GeV 


> 15 GeV 


|?7(isol fj)\ 






<«0.8 


< 1.7 


p T (tag n) 




> 4 GeV 




> 4 GeV 


|r?(tag pt)\ 




< 1.7 




< 1.7 


E T (ietl) 


> 20 GeV 


> 15 GeV 


> 25 GeV 


> 15 GeV 


E T (jet2) 


> 15 GeV 


> 10 GeV 


> 15 GeV 


> 10 GeV 


£ T (jet3) 


> 15 GeV 


> 10 GeV 


> 15 GeV 


> 5 GeV 


£ r (jet4) 






> 15 GeV 


> 5 GeV 


k(jetl)| 


< 2.5 


< 2.5 


< 3.0 


< 3.0 


k(jet2)| 


< 3.0 


< 3.0 


< 4.0 


< 4.0 


k(jet3)| 


< 3.0 


< 3.0 


< 4.0 


< 4.0 


|r ? (jet4)| 






< 4.0 


< 4.0 


No. of jets 


2 or 3 


2 or 3 


2 to 4 


2 to 4 


Et 


> 20 GeV (e in CC) 


> 15 GeV (e in CC) 


> 20 GeV 


> 15 GeV 




> 25 GeV (e in EC) 


> 20 GeV (e in EC) 






NN cosm i c 








Output > 0.7 
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The analysis starts with the simple selection criteria listed in Table |. The requirements 
are determined by comparing distributions for the summed backgrounds with those of the 
data before most of the initial selection criteria are applied, and rejecting regions where 
there is poor agreement between data and the model for the sum of signals and backgrounds. 
The backgrounds are modeled using data weighted to represent the mis-ID e or mis-ID \i 
backgrounds, and six samples of Monte Carlo (MC) events: ti, Wbb, Wcc (including Wcs 
and Wss), Wjj (where j represents u, d, or g), WW, and WZ production. 

The final requirement in Table | on the /i+jets/tag decay channel is a cutoff on the output 
of a neural network trained to reject events in which a cosmic ray has been misidentified as 



an isolated muon. The network uses the MLPfit package [13], and has seven input nodes, 
15 hidden nodes, and one output node. The input variables are described in Table |J. The 
pseudo-three-dimensional impact parameter is defined as IP3d = ylPly + IPnb where "BV" 
stands for "bend view" and "NB" for "non-bend" view for the muon trajectory through the 
spectrometer toroids. 

TABLE II. Input variables to the cosmic ray neural network. 

Variables Description 
A0(/x, tag /i) Opening angle in the transverse plane between the high-pp 

muon and tagging muon 
Pr(lJ'), PT(tag fi) Muon transverse momentum 

2 V ert(^)) -Zvcrt(tag fi) Position of the primary vertex projected from the muon's 

track in the calorimeter 
IP3d(/u), IP3d(tag fi) Pseudo-3d impact parameter of the muon's trajectory 

relative to the beam axis 



The cosmic-ray-rejection network (NN cosm i C ) is trained on a background sample of 575 
data events chosen to contain cosmic rays by requiring A</>(/z, tag [i) > 2.4, and on an equal- 
sized cosmic-ray-free signal sample of MC tb, tqb, ti, and ly+jets events. In the background 
training sample, the first muon passes either the isolated or nonisolated identification criteria. 
The results of the training are shown in Fig. |l|(a). We accept events if the value of the 
network output is greater than 0.79, a cutoff designed to maximize rejection of the cosmic-ray 
component of the background. This selects 73% of the s-channel single-top-quark acceptance, 
76% of the t-channel acceptance, 47% of the background included in the model (i.e., not 
including the rejected cosmic-ray component), and 28% of the data. The result of applying 
the network to the data and to the model for signal+background is shown in Fig. |TJ(b) . 

After applying the initial selection criteria, the combined acceptance for the s-channel 
tb signal is 3.8%; for the t-channel tqb signal it is 3.6%. The signal-to-background (S:B) 
ratios range from 1:40 to 1:470, depending on the production and decay channels. In what 
follows, we improve on the S:B ratios significantly by using arrays of neural networks that 
reject background while retaining adequate signal acceptance. The networks are trained to 
recognize detailed features of the signals and backgrounds, including correlations among the 
kinematic variables. They thereby provide superior separation of signal from background 
relative to classical selection techniques, where no correlations are taken into account. 
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Without the neural networks, no useful information about single-top- quark production can 
be obtained from the untagged events because the S:B ratio is so poor. Using the neural 
networks allows the untagged channels to provide as much sensitivity as the tagged ones. 




NNcosmic Output NN cosmic Output 



FIG. 1. Output from the neural network used to reject cosmic-ray contamination in the 
tagged muon+jets decay channel: (a) shows the results for the two training samples, and 
(b) shows the output for the sum of the signals and modeled backgrounds, and the data. 



For training the neural networks, we divide the background event samples into five sets. 
There are three sets for VF+jets events: u Wjj" u Wbb" (for the combined Wbb and Wcc MC 
samples), and lt WW" (for the combined WW and WZ MC sets); and a set each for the ti 
and misidentified-lepton ( "mis-ID /" ) backgrounds. Each analysis has five parallel networks, 
one to reject each type of background. There are four separate analyses: tb — ► e + jets, 
tb — ► /Ji + jets, tqb — > e + jets, and tqb — > \i + jets. We use the same networks for both 
untagged and tagged events combined, but choose different cutoffs on the output variables, 
depending on whether there is a tagging muon. This provides eight sets of results from 20 
neural networks and 40 cutoffs on the outputs. 



We use the package MLPfit |T3[ , which has multi-layered perceptrons with a feed- forward 
structure and back-propagation of the errors for efficient computation. The networks are 
trained on samples of 2000-9000 signal MC events and background sets of the same size, 
using the "hybrid linear Broyden-Fletcher-Goldfarb-Shanno" learning method |JL4{| . The 
performance of a neural network of course depends on the choice of input variables, which 
should be selected to provide maximal discrimination between signal and background, and 
should reflect as many different properties of the events as possible. However, too many 
input variables can worsen performance if the additional information is weak compared to 
the noise they introduce into the analysis. The variables chosen as inputs to the neural 
networks are defined in Table |TJ. 

The "best jet" in an event is the one that, when combined with the isolated lepton and 
the neutrino, generates an invariant mass closest to that of the top quark (174.3 GeV). 
The momentum components of the neutrino are derived from I^t by assuming that it 
and the lepton arise from the decay of a IF boson. Of the two possible solutions to the 
quadratic relation for the W mass, the one with the smallest absolute value of the neutrino's 
longitudinal momentum is chosen. The invariant mass variable Mb es t in Table |TXT] uses the 
best jet in its definition. Several variables use all the jets in the event ("alljets"). Those 
which exclude the best jet are denoted with a "prime". 
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TABLE III. Input variables to the neural networks. 



Variable 


Description 


Wjj 


Wbb 


WW 


tt 


Mis-ID / 


4 etl 

1 


Transverse energy of the highest- Et iet 


V 


V 


V 




v 7 

V 


i- 1 jet2 

JZ/rp 


Transverse energy of the second-highest--E/r jet 


V 


V 


V 




V 


1 r?; ot i 1 

1 /JCL1 | 


Absolute pseudorapidity of iet 1 


V 


V 






v 7 

V 


1 7 ?jct2 1 


Absolute pseudorapidity of iet 2 


V 


V 










Transverse momentum of the jetl-jet2 system 


V 


V 


V 




v 7 

V 




Scalar sum of the transverse energies of jet 1 and jet 2 






V 






1 


Transverse mass of the jetl-jet2 system 


V 


V 


v 7 

V 








Invariant mass of the jetl-jet2 system 


V 

V 


v 7 

V 






v 7 

V 


AR^ 2 


Opening angle between iet 1 and iet 2 


V 


V 








|yjij2| 


Absolute rapidity of the ietl-iet2 system 


V 

V 


V 

V 


v 7 

V 




v 7 

V 


alljets 
Pt 


Transverse momentum of the alljets system 










v 7 


-^alljets 


Invariant mass of the alljets system 








v 7 




Man 


Invariant mass of the lepton-neutrino- alljets system (= \f§) 


V 


V 




v 7 


v 7 


M bcst 


Invariant mass of the lepton-neutrino-bestjet system 


V 


V 


v 7 




v 7 


rjall 

min 


Smallest eigenvalue of the momentum tensor for the i-^-alljets system 


V 


I 

V 


v 7 




v 7 


/ r-,ietl V 

lEif 1 


Transverse energy of the highest- Et jet that is not the best jet 








v 7 






Transverse energy of the second-highest-£r jet that is not the best jet 








v 7 






Vector sum of the transverse energies of all the jets, except for the best jet 






v 7 


v 7 


v 7 


/^alljetsV 

V 1 J 


Scalar sum of the transverse energies of all the jets, except for the best jet 








v 7 




(-^alljets) 


Scalar sum of the energies of all the jets, except for the best jet 








v 7 




(-^alljets)' 


Invariant mass of all the jets, except for the best jet 








v 7 




p¥ 


Transverse momentum of the lepton-neutrino system 






v 7 








Transverse mass of the lepton-neutrino system 










v 7 


W alliets 

Pt ~Pt 


Difference between transverse momenta of the l-v and alljets systems 






v 7 




v 7 




M w -M^ 2 
M w 


Fractional difference between invariant masses of the l-v and jetl-jet2 systems 






v 7 






Et 


Missing transverse energy 








v 7 


v 7 


tag/il 

Pt 


Transverse momentum of the highest-pr tagging muon 


V 


V 


v 7 




v 7 


tag/i2 

Pt 


Transverse momentum of the second-highest-pT tagging muon 


V 


V 


v 7 




v 7 



We optimize the performance of each network by choosing the number of hidden nodes 
that minimizes the network's error function. The numbers of nodes are given in Table IV. 



TABLE IV. Number of input (i), hidden (h), and output 
(o) nodes (i-h-o) for each of the neural networks. 



Background 


s-channel tb 


t-channel tqb 


Set 


e+jets 


/i+jets 


e+jets 


/i+jets 


Wjj 


14-16-1 


14-24-1 


14-18-1 


14-21-1 


Wbb 


14-19-1 


14-19-1 


14-20-1 


14-17-1 


WW 


14-30-1 


14-19-1 


14-27-1 


14-23-1 


ti 


9-20-1 


9-20-1 


9-20-1 


9-19-1 


Mis-ID I 


16-17-1 


16-20-1 


16-28-1 


16-15-1 



Figure |2] illustrates the outputs of one set of five networks, after training on t-channel 
single-top-quark signals and each of the five background sets in the combined untagged and 
tagged electron+jets decay channels. The least amount of discrimination is obtained for 
the Wjj events, which is unfortunate since this process has a large cross section. Also, a 
significant fraction of the tqb signal cannot be differentiated from ti background. (This is 
not the case for the lower jet-multiplicity s-channel tb events.) The best separation of signal 
from background is obtained for events with a misidentified electron. 

The outputs of the neural networks for the tqb search in the untagged electron+jets decay 
channel are shown in Fig. [3|. They are obtained by passing all the signal and background 
events, and then the data, through each network. The cutoffs on the outputs are simulta- 
neously chosen by minimizing the expected limit on the cross section in this channel. 

The signal acceptances and numbers of events predicted to remain in the data after the 
initial event selection criteria are shown in Table [V|. Table |V| shows the acceptances and 
event yields after the neural network selections. We measure the signal acceptances using 
MC samples of s-channel and t-channel single-top-quark events from the compHEP event 
generator |I5|, with the pythia package ]T6| used to simulate fragmentation, initial-state 



and final-state radiation, the underlying event, and leptonic decays of the W boson. For 
all event samples, the PDF is CTEQ3M W71. The MC events are processed through a 
detector simulation program based on the geant package [HJ and a trigger simulation, and 
are then reconstructed. We apply all selection criteria directly to the reconstructed MC 
events, except for several particle identification requirements that are taken into account 
using factors obtained from other D0 data. 

We calculate the acceptance for ti pairs and for the five subprocesses for l^+jets in a 
manner similar to that used for signal, and then convert to a number of events using the 
integrated luminosity for each channel and the appropriate cross section. The ti background 



is modeled using herwig |19|1 . The Wbb, Wcc, and Wjj processes use CompHEP, followed by 



pythia, and the diboson processes are from pythia. The cross sections are DCs measured 



value for tt |20 ] , leading-order values for Wbb, Wcc, and Wjj pl[ , and next-to- leading-order 



values for WW and WZ 
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NN Output NN Output 

FIG. 2. Results of training the five neural networks used to separate t-channel tqb signal from 
background in the untagged and tagged electron+jets decay channels. For each plot, the lighter 
curves are for background and the darker curves for the signal. The solid curves are for untagged 
events, and the dashed curves for tagged ones. 




NN Output NN Output 

FIG. 3. Outputs from the five neural networks used to separate t-channel tqb signal from 
background in the untagged electron+jets decay channel. For each plot, the upper curve with 
error band shows the sum of signal and all backgrounds, the lower curve shows 50 times the 
expected single-top-quark signal, the shaded circles with error bars show the data after the initial 
event selection, and the shaded histogram shows the data after all neural network selection criteria 
have been applied, except for the network shown in that plot. 
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TABLE V. Signal acceptances (as percentages of the total cross sections) and numbers 
of events expected after application of initial selection criteria. 



e+jets/notag e+jets/tag /Lt+jets/notag /Li+jets/tag 

Signal Acceptances 

tb (2.04 ± 0.14)% (0.32 ± 0.03)% (1.27 ± 0.12)% (0.16 ± 0.01)% 

tqb (1.81 ± 0.15)% (0.20 ± 0.02)% (1.44 ± 0.14)% (0.12 ± 0.01)% 

Numbers of Events 



tb 


1.41 ± 0.25 


0.22 ± 0.04 


0.84 ± 0.16 


0.11 ± 0.02 


tqb 


2.44 ± 0.43 


0.27 ± 0.05 


1.85 ± 0.34 


0.16 ± 0.03 


tt 


6 ± 2 


1.42 ± 0.43 


8 ± 3 


1.45 ± 0.44 


Wbb 


4 ± 1 


0.69 ± 0.23 


2 ± 1 


0.27 ± 0.09 


Wcc 


17 ± 6 


0.80 ± 0.27 


11 ± 4 


0.24 ± 0.09 


Wjj 


499 ± 120 


2.88 ± 0.95 


332 ± 84 


0.65 ± 0.28 


WW 


17 ± 2 


0.31 ± 0.05 


13 ± 2 


0.24 ± 0.05 


wz 


3 ± 


0.11 ± 0.02 


2 ± 


0.06 ± 0.01 


Mis-ID I 


112 ± 13 


9.79 ± 1.32 


19 ± 4 


2.06 ± 1.02 


Bkgd for tb 


661 ± 130 


16.28 ± 2.03 


389 ± 91 


5.93 ± 1.22 


Bkgd for tqb 


660 ± 130 


16.23 ± 2.02 


388 ± 91 


5.88 ± 1.22 


Data 


558 


14 


398 
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In the electron channel, the misidentified-electron background is measured using multijet 
data. For each jet that passes the electron E? and \r)\ requirements, the events are weighted 
by the probability that a jet mimics an electron. These probabilities are determined from 
the same multijet sample, but for Et < 15 GeV, and are found to be (0.0231 ± 0.0039)% 
(CC), and (0.0850 ± 0.0118)% (EC) in untagged events, and (0.0154 ± 0.0019)% (CC), and 
(0.0612 ± 0.0057)% (EC) in tagged events. The probabilities are independent of the jet E T 
for Et > 20 GeV. We normalize the integrated luminosity of the multijet sample to match 
the data sample used in the search for signal, and correct for a small difference in trigger 
efficiency between the two samples. 

In the muon channel, the misidentified-muon background is from bb events; it is produced 
when one or both of the b quarks decays semileptonically to a muon, and one muon is 
misidentified as isolated. There are two ways such events can mimic signals. First, one 
of the b jets may not be reconstructed, and its muon can therefore appear to be isolated. 
Second, a muon can be emitted wide of its jet and be reconstructed as an isolated muon. 
The background from each source is measured using data collected with the same triggers as 
used for the muon signal. The events are required to pass all selection criteria, except that 
the muon, which otherwise passes the isolated muon requirements, is within a jet. Events 
with truly isolated muons are excluded. Each event is then weighted by the probability that 
a nonisolated muon is reconstructed as an isolated one. This probability is measured using 
the same data, but for Et < 15 GeV, and is found to be a few percent for each source on 
average. The probabilities are parametrized as a function of the muon pt] they are higher at 
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low pt and fall to zero for px > 28 GeV. We calculate a weighted average of the two results 
to obtain the number of mis-ID /i background events. 

TABLE VI. Signal acceptances (as percentages of the total cross sections) and numbers 
of events expected after application of both the initial and neural network selection criteria. 



e+jets/notag e+jets/tag /i+jets/notag /i+jets/tag 



s-channel tb Search 








Signal Acceptance 






± 




tb 


(0.29 


± 


0.02)% 


(0.20 


± 0.02)% 


(0.24 ± 


0.03)% 


(0.12 


± 


0.01)% 












Numbers of Events 










tb 


0.20 


± 


0.04 


0.14 


± 0.03 


0.16 ± 


0.03 


0.08 


± 


0.02 


tqb 


0.15 


± 


0.03 


0.10 


± 0.02 


0.14 ± 


0.03 


0.06 


± 


0.01 


ti 


0.29 


± 


0.10 


0.20 


± 0.06 


0.26 ± 


0.08 


0.10 


± 


0.03 


Wbb 


0.15 


± 


0.05 


0.30 


± 0.10 


0.16 ± 


0.06 


0.16 


± 


0.06 


Wcc 


0.77 


1 

± 


0.27 


0.32 


± 0.11 


0.74 ± 


0.27 


n i o 

0.12 


± 


0.05 


Wjj 


13.76 


± 


3.45 


0.89 


± 0.42 


14.56 ± 


4.28 


0.43 


± 


0.23 


WW 


0.24 


± 


0.11 


0.08 


± 0.03 


0.19 ± 


0.07 


0.12 


± 


0.03 


wz 


0.06 


1 


0.02 


0.04 


± 0.01 


0.06 ± 


r\ no 

0.02 


n no 

0.02 


1 


0.01 


Mis-ID I 


1.06 


± 


0.23 


0.37 


± 0.11 


0.00 ± 


0.00 


n nn 

0.06 


1 

± 


0.05 


Total Bkgd 


16.49 


± 


3.83 


2.29 


± 0.61 


16.10 ± 


4.66 


1.07 


± 


0.32 


Data 
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t-channel tqb Search 








Signal Acceptance 






± 




tqb 


(0.29 


± 


0.03)% 


(0.13 


± 0.01)% 


(0.38 ± 


0.04)% 


(0.08 


± 


0.01)% 












Numbers of Events 










tqb 


0.38 


± 


0.08 


0.17 


± 0.03 


0.50 ± 


0.10 


0.11 


± 


0.02 


tb 


0.10 


± 


0.02 


0.12 


± 0.02 


0.10 ± 


0.02 


0.07 


± 


0.01 


ti 


0.78 


± 


0.24 


0.46 


± 0.14 


1.63 ± 


0.50 


0.20 


± 


0.06 


Wbb 


0.07 


± 


0.03 


0.20 


± 0.07 


0.08 ± 


0.03 


0.11 


± 


0.04 


Wcc 


0.36 


± 


0.13 


0.23 


± 0.09 


0.42 ± 


0.16 


0.10 


± 


0.04 


Wjj 


10.16 


± 


3.26 


0.83 


± 0.39 


13.18 ± 


4.60 


0.22 


± 


0.13 


WW 


0.43 


± 


0.11 


0.04 


± 0.01 


1.08 ± 


0.31 


0.09 


± 


0.03 


WZ 


0.08 


± 


0.03 


0.03 


± 0.01 


0.20 ± 


0.05 


0.03 


± 


0.01 


Mis-ID I 


0.76 


± 


0.20 


0.29 


± 0.08 


0.03 ± 


0.03 


0.08 


± 


0.03 


Total Bkgd 


12.75 


± 


3.58 


2.22 


± 0.56 


16.73 ± 


5.13 


0.91 


± 


0.23 


Data 
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After applying the neural network selection criteria, the combined acceptances are 0.86% 
for the tb signal and 0.88% for the tqb signal, with the S:B ratios increased to between 1:9 
and 1:99, a factor of 4-8 improvement compared with those after the initial event selections. 
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We use a Bayesian approach p4| to calculate limits on the cross sections for single-top- 



quark production in the s-channel and t-channel modes. The inputs are the numbers of 
observed events, the signal acceptances and backgrounds, and the integrated luminosities. 
Covariance matrices are used to describe the correlated uncertainties on these quantities. A 
flat prior is used for the single-top-quark cross section, and a multivariate Gaussian prior 
for the other quantities. We calculate the likelihood functions in each decay channel and 
combine them to obtain the following 95% confidence level upper limits: 

• a(pp -> tb + X) < 17 pb 

• a(pp -> tqb + X) < 22 pb. 

The contributions of each decay channel to these results are shown in Table |V11| . 



TABLE VII. The 95% confidence level upper limits on cross sections for the two 



production modes of 


single top quarks. 


Values are 


in picobarns. 






Initial Selections 


Neural Network Selections 




e+jets 


//+jets 


e+jets 


e+jets 


/x+jets 


e+jets 


Channel 






+ //+jets 






+ /x+jets 


s-channel tb 














untagged 


118 


165 


107 


44 


45 


35 


tagged 


27 


108 


35 


26 


39 


19 


untagged+tagged 


27 


104 


36 


22 


26 


17 


t-channel tqb 














untagged 


131 


156 


141 


41 


43 


33 


tagged 


43 


141 


55 


43 


59 


30 


untagged+tagged 


42 


128 


60 


27 


32 
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To conclude, we have searched for the electroweak production of single top quarks using 
a neural-network signal-selection technique. We find no evidence for such production and 
set upper limits on the cross sections for s-channel production of tb and t-channel production 
of tqb. The limits are consistent with expectations from the standard model. 
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single-top-quark cross sections for us. We thank the staffs at Fermilab and collaborating 
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Departments of Atomic Energy and Science and Education (India), Colciencias (Colombia), 
CONACyT (Mexico), Ministry of Education and KOSEF (Korea), CONICET and UBACyT 
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